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Abstract. We consider an operator A 2 + A{x) ■ D + q(x) with the Navier 
boundary conditions on a bounded domain in M. n , n > 3. We show that a first 
order perturbation A(x) ■ D + q can be determined uniquely by measuring the 
Dirichlet-to-Neumann map on possibly very small subsets of the boundary of 
the domain. Notice that the corresponding result does not hold in general for 
a first order perturbation of the Laplacian. 



1. Introduction and statement of results 

Let Q C M. n , n > 3, be a bounded simply connected domain with C°° connected 
boundary, and let us consider the following equation, 



£A q {x, D)u — in Q, 



where 



Ca, q {x, D) = A 2 + ^2 Aj(x)D,j + q(x) = A 2 + A(x) ■ D + q(x), 

3=1 

D = r x V, A = {Aj)i<j< n G C 4 (n,C n ), and q G L°°(fi,C). The operator C A , g , 
equipped with the domain 

V{C A>q ) = {ue H 4 (n) : u\ dn = (Au)\en = 0}, (1.1) 

is an unbounded closed operator on L 2 (Q) with purely discrete spectrum, see 
[14] . The boundary conditions in (11.11) are called the Navier boundary conditions. 
Physically, the operator C A , q with the domain T>(C A i9 ) arises when considering the 
equilibrium configuration of an elastic plate which is hinged along the boundary, 
see [S]. Let us make the following assumption, 



(A) is not an eigenvalue of C Atq (x, D) : V(CA, q ) — > L 2 (Q). 

l 
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Under the assumption (A), for any (/ , fx) G H 7 ^ 2 (dQ) x H 3 ^ 2 (dQ), the boundary 
value problem 

C-A,q u = in Q, 

u = f on dQ, (1.2) 
Am = fx on dQ, 

has a unique solution u G if 4 (f2). Let z/ be the unit outer normal to the boundary 
dQ. We then define the Dirichlet-to-Neumann map H A} q by 

M Aa ■ H 7/2 (dQ) x H 3/2 (dQ) -+H 5/2 (dQ) x H 1/2 (dQ), 
J^A, q (fo,fi) =(d u u\ dn ,d u (Au)\ dn ) 

where u G H 4 (Q) is the solution to the problem fll.2p . Let us also introduce the 
set of the Cauchy data C AtQ for the operator C A ^ q defined as follows, 

C A , q = {{u\dn, {Au)\ 9 n,d v u\an,d v (Au)\ d n) : u G H 4 (Q),£ A>q u = in O}. 

When the assumption (A) holds, the set C A>q is the graph of the Dirichlet-to- 
Neumann map M A , q . 

It was shown in [23] that the first order perturbation A(x) ■ D + q(x) of the 
polyharmonic operator (— A) m , m > 2, can be recovered from the knowledge 
of the set of the Cauchy data C AtQ on the boundary of Q. Notice that for the 
corresponding problem for a first order perturbation of the Laplacian, this is no 
longer true, due to the gauge invariance of boundary measurements, see [T9J [30j 
[37] . In this case, the first order perturbation can be recovered only modulo a 
gauge transformation, [30| [37]. 

In this paper we are concerned with the inverse problem of determining the first 
order perturbation A(x)-D + q(x) of the biharmonic operator from the knowledge 
of the Dirichlet-to-Neumann map Af A>q , given only on a part of the boundary 
dQ. 

In many applications, say, arising in geophysics, performing measurements on the 
entire boundary could be either impossible or too cost consuming. One is there- 
fore naturally led to an inverse boundary value problem with partial measure- 
ments. Substantial progress has been made recently on the partial data problems 
in the context of electrical impedance tomography as well as the Schrodinger 
equation. Specifically, in the paper [3] it was shown that an unknown conductiv- 
ity is determined uniquely by performing voltage-to-current measurements on, 
roughly speaking, a half of the boundary. The main technical tool was a boundary 
Carleman estimate with a linear weight. The Carleman estimates approach to the 
partial data problem was very much advanced in the work [20] • Here, rather than 
working with linear weights, a broader class of limiting Carleman weights was in- 
troduced and employed. The work [20J contains a global uniqueness result for the 
conductivity equation, assuming that the voltage-to-current map is measured on 
a possibly very small subset of the boundary, with the precise shape depending 
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on the geometry of the boundary. The limiting Carleman weights approach of 
[2Ti] has led to subsequent important developments in the partial data problem 
for the magnetic Schrodinger operator [7j, [21] and the Dirac system |36j. To 
the best of our knowledge, applications of this approach to partial data problems 
for other important equations and systems of mathematical physics have not yet 
been explored. The purpose of this paper is to apply the techniques of Carleman 
estimates to the partial data problem for the perturbed biharmonic operator. 

We should also mention another approach to the partial data problems for the 
conductivity equation, which is due to [18], and which is based on reflection 
arguments. In this approach, the subset of the boundary, where the measurements 
are performed is such that the inaccessible part of the boundary is a subset of a 
hyperplane or a sphere. The work [5] gives a partial data result analogous to [TS] 
for the Maxwell equations. 

Let us now proceed to describe the precise assumptions and results. Let x G 
]R n \ ch(Q), where ch(f2) is the convex hull of Q. Following [20], we define the 
front face of dQ with respect to xq by 

F(x ) = {x G dQ : (x - x ) ■ v(x) < 0}, (1.3) 

and let F be an open neighborhood of F(xo) in dQ. The main result of this paper 
is as follows. 

Theorem 1.1. Let Q C M. n , n > 3, be a bounded simply connected domain with 
C°° connected boundary, and letA^,A^ G C 4 (n,C n ) andq^\q^ G L°°(Q,C), 
be such that the assumption (A) is satisfied for both operators. If 

MAwMfo>h)\F=MAwMfoJi)\F f° r aU (A, A) e H 7 l\dQ) x HV 2 (dQ), 
then = A^ and q {l) = q {2) in Q. 

Following [20] , we say that an open set Q C M n with smooth boundary is strongly 
star shaped with respect to a point Xi G dQ, if every line through x\ which 
is not contained in the tangent hyperplane cuts the boundary dQ at precisely 
two distinct points X\ and X2, with transversal intersection at x<i- We have the 
following corollary of Theorem 11.11 

Corollary 1.2. Let Q C MJ 1 , n > 3, be a bounded simply connected domain with 
C°° connected boundary, and let x\ G dQ be such that the tangent hyperplane of 
dQ at X\ only intersects dQ at x\ and Q is strongly star shaped with respect to 
x x . Furthermore, let A^,A^ G C 4 (n,C n ) and qW,qW G L°°(Q,C), be such 
that the assumption (A) is satisfied for both operators. If for a neighborhood F 
of Xi in dQ, we have 

■^AW, q W (fo, A) | j? = NaWjW (fo, fi) I j? for all (f , A) e H 7 ' 2 {dQ) x H^ 2 {dQ), 
then A^ = A^ and = q {2) in Q. 



4 



KRUPCHYK, LASSAS, AND UHLMANN 



Notice that if Q is strictly convex, then the assumptions on Q of Corollary II .21 are 
satisfied for any x\ £ dfl, and therefore, measuring the Dirichlet-to-Neumann 
map on an arbitrarily small open subset of the boundary determines the first 
order perturbation uniquely. 

Let F be an open neighborhood of the front face of dQ with respect to Xq, defined 
in ( 11. 3j) . Associated to F, consider the set of the Cauchy data for the first 
order perturbation of the biharmonic operator, which is based on the Dirichlet 
boundary conditions (u\gn, d v u\Qn), 

C Ati), q U) = {( u \an,duu\dn,dlu\ dn ,dlu\p) : u £ # 4 (0), C A u), q U)U = in Q}, 

j = 1,2. Notice that the Dirichlet boundary conditions correspond to the clamped 
plate equation. We have the following partial data result. 

Corollary 1.3. Let Q C M. n , n > 3, be a bounded simply connected domain with 
C°°jconnected [boundary, and letA^,A^ £ C 4 (fi,C n ) and , g (2) £ L°°(fi,C). 

V c lw, q w = c lv), g w> then A(1) = A(2) and = <l (2) m Q - 

Corollary 1 1 . 31 follows from Theorem 1 1 . 1 1 and the explicit description for the Lapla- 
cian in the boundary normal coordinates, see [26J. 

Finally, let us mention that the study of inverse boundary value problems has 
a long and distinguished tradition, in particular, in the context of electrical 
impedance tomography, see [Tj [2j [29J [38] for the two dimensional case, and 
[5], Dm EEJ HS1 [39] for the case of higher dimensions, as well as in inverse bound- 
ary value problems and inverse scattering problems for the Schrodinger equation 
[31 El EE El EDI ED EH EZ] , and in elliptic inverse problems on Riemannian man- 
ifolds, [151 [161 HU [211 [251 [26] . F° r sufficiently non-regular coefficients, the inverse 
problems are not uniquely solvable, see [101 [TTJ , even when the measurements are 
performed on the whole boundary. These counterexamples are closely related to 
the so-called invisibility cloaking, see e.g. [TJl [131 [221 [23 132] . 

The paper is organized as follows. In Section 2 we construct complex geometric 
optics solutions for the perturbed biharmonic operator, using the methods of 
Carleman estimates with limiting Carleman weights. Section 3 is devoted to 
Carleman estimates with boundary terms for the perturbed biharmonic operator. 
These estimates are crucial when estimating away the boundary terms in the 
derivation of the main integral identity, which is carried out in Section 4. The 
final Section 5 is concerned with the determination of the first order perturbation, 
relying upon the main integral identity. We notice that the unique identifiability 
of the vector field part of the perturbation becomes possible thanks to special 
properties of the amplitudes in the complex geometric optics solutions. 
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2. Construction of complex geometric optics solutions 

Let Q C R™, n > 3, be a bounded domain with C°°-boundary. Following [7J 120] . 
we shall use the method of Carleman estimates to construct complex geometric 
optics solutions for the equation Ca, q u = in Q, with A G C 4 (fi,C n ) and q G 
L°°(tt,C). 

First we shall derive a Carleman estimate for the semiclassical biharmonic oper- 
ator (— h 2 A) 2 , where h > is a small parameter, by iterating the corresponding 
Carleman estimate for the semiclassical Laplacian —h 2 A, which we now proceed 
to recall following [2"0| 133] . Let Q be an open set in R n such that Q CC Q and 
ip G C°°(f2,R). Consider the conjugated operator 

P v = e %(-h 2 A)e-% 

and its semiclassical principal symbol 

P v (x,Z) = ¥ + - |V<^| 2 , i6fi, ^Gl". (2.1) 

Following [20], we say that tp is a limiting Carleman weight for —h 2 A in f2, if 
V<£> 7^ in Q and the Poisson bracket of Rep v and lmp v satisfies, 

{Rep^jImp^Xx,^) = when p<p(x,f)=0, (x,0efixlR n . 

Examples are linear weights </?(£) = a-x, a G R™, |a| = 1, and logarithmic weights 
ip(x) = log | a; — xq\, with xq Q. In this paper we shall use the logarithmic 
weights. 

In what follows we shall equip the standard Sobolev space H s (M. n ), s G R, with 
the semiclassical norm ||w||# s s cl = || (hD) s u\\L2. Here (£) = (1 + l^] 2 ) 1 / 2 . We shall 
need the following result, obtained in [35], generalizing the Carleman estimate 
with a gain of one derivative, proven in [20] . 

Proposition 2.1. Let if be a limiting Carleman weight for the semiclassical 
Laplacian on Q. Then the Carleman estimate 

||e*(-/i 2 A)e-*u||^ > -^-|M|^ + 2, C SiQ >0, (2.2) 

C s ,n scl 

holds for all u G C£°(Q), s G R and all h > small enough. 

Iterating the Carleman estimate (12.21) two times, we get the following Carleman 
estimate for the biharmonic operator, 

\\e%(h 2 A) 2 e-Zu\\ H!ci > J^\\ u \\ H s +ij (2.3) 

C s ,n scl 

for all u G C£°(Q), s G R and h > small. 
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Let A E W 1 ' 00 ^,^ 1 ) and q E L°°(Q, C). Then to add the perturbation /i 4 g to 
the estimate (I2.3p . we assume that —4 < s < and use that 

scl sc \ 

To add the perturbation 

h 3 e%(A ■ hD)e~% = h 3 (A -hD + iA- V<p) 
to the estimate f)2.3p . assuming that —4 < s < 0, we need the following estimates 
\\{A- V<p)u\\ H!ci < \\A-V(p\\ L °°\\u\\ H s+4, 

n 

\\A ■ hDu\\ HLi < ^ WhD^Aju)^ + 0(h)\\(divA)u\\ H s ci 

3=1 

n 

< 0(1) WAiuWjr+i + 0(h)\\u\\ H s V < 0(1)|| W ||^ 4 . 

' * scl scl scl 

3=1 

When obtaining the last inequality, we notice that the operator, given by mul- 
tiplication by Aj, maps H^ 4 — > . To see this by complex interpolation it 
suffices to consider the cases s = and s = —4. 

Let 

C ip = e^h L A , q e h . 

Thus, we obtain the following Carleman estimate for a first order perturbation 
of the biharmonic operator. 

Proposition 2.2. Let A E W 1,0 °(fi, C n ), q E L°°(tt,C), and <p be a limiting 
Carleman weight for the semiclassical Laplacian on Q. If —4 < s < 0, then for 
h > small enough, one has 

h 2 

\\C v u\\ H * cl > IMIh»+*> (2.4) 

for alluE C °°(fi). 

The formal L 2 -adjoint of C v is given by 

= e~^(h 4 C Ati -i v . A+ g)eh. 

Notice that if <p is a limiting Carleman weight, then so is — (p. This implies that 
the Carleman estimate f!2.4[) holds also for the formal adjoint C* . 

To construct complex geometric optics solution we need the following solvability 
result, similar to [7j. The proof is essentially well-known, and is included here for 
the convenience of the reader. 
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Proposition 2.3. Let A G W lt00 (Q, C n ), q G L°°(f2,C) ; and let be a limiting 
Carleman weight for the semiclassical Laplacian onQ. If h > is small enough, 
then for any v G L 2 (Q), there is a solution u G H 4 (Q) of the equation 

C v u = v in Q, 

which satisfies 

„ „ c„ „ 

\\u\\ HL < J^\\V\\ L 2. 

Proof. Consider the following complex linear functional 

l : c;c™(n) -> c, c; w ^(w,v) L2 . 

By the Carleman estimate (12. 4p for the formal adjoint £*, the map L is well- 
defined. 



Let w G C7 °°(fi). We have 

1 -JMrol <C II in\\ r n II 1 ) II t-1 <T 

h- " ' ""-i 



c 



showing that L is bounded in the if _4 -norm. Thus, by the Hahn-Banach theorem, 
we may extend L to a linear continuous functional L on H~ A (R n ) without increas- 
ing the norm. By the Riesz representation theorem, there exists u G H 4 (M. n ) such 
that for all w G H~ 4 (R n ), 

~ C 
L(w) = (w,tt)(H~*,ff4), and \\u\\ H ^ < j»\\v\\ip. 

Here (•, -)(h- 4 ,h 4 ) stands for the usual L 2 -duality. It follows that C v u = v in Q. 
This completes the proof. □ 

Our next goal is to construct complex geometric optics solutions of the equation 

£-A, q u = in Q, (2.5) 
with A G C 4 (f2, C 4 ) and q G L°°(f2, C), i.e. solutions of the following form, 

u(x; h) = e h (a (x) + ha\(x) + r(x; h)). (2.6) 

Here (p G C°°(f2,R) is a limiting Carleman weight for the semiclassical Laplacian 
on O, ip G C°°(Q, R) is a solution to the eikonal equation p^(x, Vip) =0 in Q, 
where p v is given by (12.11) . i.e. 

|V?/f = |Vyf, V<^-V^ = 0, in n, (2.7) 

the amplitudes a G C°°(fi) and a\ G C 4 (f2) are solutions of the first and second 
transport equations, and r is a correction term, satisfying ||r|| H 4 ^ = 0{h 2 ). 
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Following [3 [20], we fix a point xq G M n \ ch(f2) and let the limiting Carleman 
weight be 

p(x) = -\og\x - x \ 2 , (2.8) 

and 

if \ n + uj-(x-xq) f X-Xq \ 
TO = o ~ arctan 77 xo / / = dls V-i I \> u , ( 2 - 9 ) 

2 y/{x - X ) 2 - {U ■ {X - X )) 2 \\X-X Q \ J 

where u G S n_1 is chosen so that if) is smooth near Q. Thus, given ip, the function 
if) satisfies the eikonal equation (12.71) near Q. 

Consider the conjugated operator, 

e^^A,^^ 4 = {h 2 A + 2hT) 2 + h 3 A-hD + h 3 A-(D<p + tDif)) + h 4 q, (2.10) 
where 

T= (V<p + iVip) ■ V + -{Ap + iAif)). (2.11) 

Substituting (I2.6P into (12.51) and collecting powers of h, we get the first transport 
equation, 

T 2 a = in ft, (2.12) 

the second transport equation 

T 2 a x = -^(AoT + To A)a - ■ {D<p + iDip)a in Q, (2.13) 

and 

e h h LA,q{e h r) = — e h ft >c^ jg (e h (a + ftaij) 

= - /i 4 A 2 (a + hat) ~ 2/> 4 (A oT + To A)a x - /i 4 A • Da 

- /i 5 A • Da! - /i 4 A • (Dip + iDif))a t - h A q(a Q + ftai). 

(2.14) 

Let us now discuss the solvability of the equations ( 12.121) . ( 12. 13ft and (12.141) . 
To this end we follow [7] and choose coordinates in R n so that Xq = and 

f2 C {x n > 0}. We set u> = e\ = (1,0k«-i), an d introduce also the cylindrical 
coordinates (xi,r9) on M. n with r > and # G § n ~ 2 . Consider the change of 
coordinates x H- (z,9), where z = x\ + ir is a complex variable. Notice that 
Imz > near Q. Then in these coordinates, we have 

7T Rez 

9? = log \ z\ = Re log 2, if) — arctan = Im log z, 

2 Imz 

when Imz > 0. Hence, 

p + iip = log 2. 

Moreover, we have 

V((p + iip) = -(ex + ie r ), (2.15) 
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where e r = (0,6), 6 E S n ~ 2 , and 

2 2(n — 2) 

V(^ + iV)-V = -^, + #) = —) =(. (2.16) 

2 z(z — z) 

In the cylindrical coordinates, the operator T, defined in (12. lip , has the form, 



2(z -f; 

Thus, it follows that (I2.12p has the form, 

9z~ ^~ 2 -l l ao = in fi. 



2(z-s) 

In particular, one can take ao E C°°(fi) satisfying ^<9 g — 2 ^"~^ ^ a = 0. The 

general solution of the latter equation is given by ao = (z — zY 2 ~ n ^ 2 go with 
g E C7°°(fi) satisfying d g g = 0. 

In the cylindrical coordinates, the second transport equation (I2.13P has the form, 

(n - 2) x 2 



d- z - ± ai = / in n, (2.17) 

2{z-z)J 

where / is given. Notice that a\ will have in general the same regularity as 
/, which is the same as the regularity of A. It follows from (I2.14p that we 
need four derivatives of a\, which explains our regularity assumptions on A, i.e. 
AeC 4 {Q,C n )- 

In order to solve (I2.17p . given / E C 4 (f2), one can find v E C 4 (f2), which satisfies 

<n-2) \ 

. 2(z-z)J { ' 

and then solve 

' (n-2). . 

°s ~ TTt =\ ) a i= v 111 

2(z — z) 

We look for a solution of (12. 18ft in the form v = e 9 vo with g E C°°(Q) satisfying 
d s g = 2 (^z) • Thus, v E C 4 (Q) can be obtained by solving d 2 v = e~ 9 f, applying 
the Cauchy transform, i.e. 

v (z,0) = — / rfReC«IniC, 

Jc C 

where \ E Cq^M" - ) is such that x = \ near Q. Hence, the second transport 
equation (12.131) is solvable globally near f2 with a solution a\ E C 4 (f2). 

Having chosen the amplitudes ao E C°°(Q) and a x E C A (Q), we obtain from 
(gUl) that 

e^h 4 C A ^(e^r) = 0(h A ), 
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in L 2 (Q). Thanks to Proposition I2.3[ for h > small enough, there exists a 
solution r G H*(Q) of (IzA such that ||r|| H 4 c] = 0{h 2 ). 

Summing up, we have the following result. 

Proposition 2.4. Let A G C A (U,C n ), q e £°°(fi,C). T/jen /or all h > small 
enough, there exist solutions u(x; h) G if 4 (fi) to £/ie equation Ca, q u = in Q, of 
the form 

u(x; h) = e >>■ (a (x) + hai(x) + r(x; h)), 

where ip G C°°(Q, R) is a limiting Carleman weight for the semiclassical Laplacian 
on tt, chosen as in (J2IBJ, ip G C°°(fi,M) is given by (|Z^)I , a G C°°(Tl) and 
ai G C 4 (fi) are solutions of the first and second transport equations (12.121) and 
(I2.13p . respectively, and r is a correction term, satisfying \\r\\ H 4 ^ = 0{h 2 ). 



3. Boundary Carleman estimates 

Following [7J [20] , in order to prove that some boundary integrals tend to zero as 
h — > 0, we shall use Carleman estimates, involving the boundary terms. 

Let Q C M n , n > 3, be a bounded domain with C°°-boundary, and Q C M n be an 
open set such that Q CC £1, and ip G C°°(f2, R) be a limiting Carleman weight 
for the semiclassical Laplacian. We define 

dtt± = {x G dtt : ±<9^(x) > 0}. 

We shall need the following result from |2Uj . 

Proposition 3.1. Let ip be a limiting Carleman weight for the semiclassical 
Laplacian on Q. Then there exists C > such that the following Carleman 
estimate 

||e"^(-/i 2 A)w|| L2 + 2 \\y/-d u (p e~hd v u\\ L 2r m ) 
> —{h\\e h u\\ H i + h 3/2 \\^/d^pe *d u u\\ L 2 {m+) ), 

holds for all u G H 2 (Q), u\qq = 0, and all h > small enough. 

Iterating (13. ip two times and adding a first order perturbation, we get the follow- 
ing boundary Carleman estimate for the first order perturbation of the biharmonic 
operator. 
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Proposition 3.2. Let A G W 1 '° C (Q, C n ), q G L°°(fi,C), and y be a limiting Car- 
leman weight for the semiclassical Laplacian on Q. Then the following estimate 

\\e^(h A C A , q )u\\ L2 + h 3 / 2 \\ y/^e-Sdvi-tfAu)]]^^ 

I _ (3.2) 

> -^{h 2 \\e-^u\\ H i + h 3/2 \\^/d^e-%d u {-h 2 Au)\\ L 2 {dn+) 

+ h 5/2 \\y/d^pe~^d v u\\ L 2 im+) ), 
holds, for all u G H 4 (Q), u\gn = {Au)\en — 0, and all h > small enough. 

Here we use that 

\\e-%h*A- Du\\ L 2 < h 3 J2 (\\ hD j( A j e ~* u )\\L* + ll(^Dj(^j e ~")HU 2 ) 

j= i V / 

< h 3 ^ ( 2 \\( hD jAj)e~%u\\ L 2 + \\AjhDj{e-%u)\\ L i + WA^D^e'^uW^ j 
j= i V / 

<0{h?)\\e-U\\ Hlci . 

Notice that when (p is given by (12. 8ft . we have d v (p{x) = ^^p^ an< ^ therefore, 
dVL_ = F(x Q ). 



4. Integral identity, needed to recover the first order 

perturbation 



We shall need the following Green's formula, see 

/ (CA, q u)vdx — / uC* Aq vdx = —i / v(x) ■ uAvdS — / d u (—Au)vdS 
Jn Jn Jan Jan 



+ / {-Au)d v vdS- / d v u{-Av)dS+ / u(d v (-Av))dS, 
'an Jan Jan 



(4.1) 



valid for all u, v G H 4 (Q). Here C* A = C-A^v-A+q is the adjoint of £-A,q, v is the 
unit outer normal to the boundary dQ, and dS is the surface measure on dfl. 

Let (/o, /i) e H 7 ' 2 (dtt) x H'^ 2 (dn) and ^ G # 4 (ft) be such that 
£>AU),qU)Uj = in fi, j = 1,2, = i^lan = /o, (Awi)|an = (Aw 2 )|an = fx- 



Then by assumption of Theorem there exists an open neighborhood F of 
F(xq) = dfl- in dfl, such that 

d„ux\p = d u u 2 \p, d„(Aux)\p = d u (Au 2 )\ F . 
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We have 

C AW , qW (u 1 -u 2 ) = (A^-A^)-Du 2 + (q^-q^)u 2 in Q. (4.2) 
Let v G H 4 (Q) satisfy C* AW q{1) v = in Q. Using (14. ip . we get 

((,4(2) _ . Du 2 )vdx + I (g {2) - q {1) )u 2 vdx 

Jn (4.3) 



_d v (-A(ux -u 2 ))vdS- / _d v {u x -u 2 )(-Av)dS. 
an\F Jan\F 

To show the equalities = A^ and = q^ 2 \ the idea is to use the identity 
(I4.3P with u 2 and v being complex geometric optics solutions and to use the 
boundary Carleman estimate (13. 2 p to show that the boundary integrals in (I4.3p . 
multiplied by some power of h, tend to zero as h — > 0. 

To construct the appropriate complex geometric optics solutions, let (p and ip be 
defined by (12. 8p and (12.91) . respectively. Then thanks to Proposition I2.4[ we can 
take 

u 2 (x- h) = e^(4 2) (x) + haf\x) + r^ 2 \x; h)), (4.4) 
v (x; h) = e^ /i (4 1) (x) + haf\x) + r«(x; h)), (4.5) 
where a { j) G C°°{Tl) and af ] G C 4 (H), j = 1, 2, are such that 



-Vyj + iVip) ■ V + ^— < J = 



Ay+jA^y (i) 



(4.6) 



and 

\\r^\\ HL = 0(h 2 ). (4.7) 
Substituting u 2 and v, given by (14. 4p and (14. 5p . in (14. 3p . we get 

(A (2) - A (1) ) • ^(Dtp + iD^)(af + ha? + r^)(a? + ^ + TW)dx 



+ / (A< 2 > - A«) ■ (Da? + /iDaS 2) + Dr^)(a? + to? 3 + r^)dx 

r _ _ _ (4. 

+ / (q (2) - g (1) )(4 2) + ^ + r^)(a? + ha? + rW)dx 



[ _d t ,(-A(u 1 -u 2 ))vdS- I _d u ( Ul -u 2 )(-Av)dS. 



'dtt\F JdU\F 

Let us now show that 
h 



/ _d„(-A(u 1 -u 2 ))vdS ->0, as h ->• +0, (4.9) 
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and 

h I d v {u x -u 2 ){-Av)dS ->0, as h -> +0, (4.10) 

where «2 and u, given by (I4.4p and (14. 5p . To this end, notice that (14. 7p implies 
that r^' = 0(1), j = 1,2, in the standard (h = 1) if 2 (fi)-norm. Hence, 

r^\ 9n = 0(l) in L\dQ), (4.11) 

and 

Vr% n = 0(1) in L 2 (dfi). (4.12) 

Moreover, it follows from (14.71) that = 0(1/ h), j = 1,2, in the standard 
H 3 (Q)-norm. Thus, 

Ar U) \ dn = 0{l/h) in L 2 {dQ). (4.13) 
Furthermore, by the definition of F(xq) and F, there exists e > fixed such that 
dtt_ = F(x ) C F £ := {xedVl: d„cp <e}cF. 

Using the Cauchy-Schwarz inequality and (14.111) . we have 

h [ d v {-A{ Ul - u 2 ))vdS 
Jan\F 

< h I \d v {-A{u x - u 2 ))\e^\a ( o ) + ha? + r {1) \dS 

Jan\F E 

<0(h)( / e|^(-A( Ml - M2 ))| 2 e^^ \\a? + ha?+r^\\ L 2 m) 
V Jan\F e J 

< 0(h)\\y^e^d u (-A (ui - w 2 ))||i2(an + ). 

By the boundary Carleman estimate (13.21) and (14.21) . we get 

0(h)\\y/d^e^d u {-A{ Ul -u 2 ))\\ L 2 iaQ+) < O{0 2 )\\e^C AWAW { Ul -u 2 )\\ L , {n) 
= O(0 2 )\\e^(AW - A«) ■ Du 2 + e^(qV - q^)u 2 \\ LHn) 

< O{0 2 )\\{A^ - A«) ■ ^{Dip + iD^af + ha? + r^)\\ L , {n) 

+ O{0 2 )\\{A^ - A«) • (L>4 2) + W?a?> + Pr( 2 ))|U 2(n) 
+ 0(^ 2 )||(g( 2 ) - g«)(4 2) + ha? + r( 2 ))|| i2(Q) < 0{h}' 2 ). 

(4.14) 

Thus, gU) follows. 
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To establish (I4.10p . we first notice that thanks to (12 .7p . we have 
Av =e^ r -^ + ^ (Q (D + ha M + r (D) 

+ 2 ' V ^ + ^ . (Va « + Waf» + Vr«) + Aa« + hAa? + Ar«.) 

This together with ( 14TH . f02j) and ( 14131) imply that 

Av| an = e^v, tJ=£>(1//i) in L 2 (dn). 
This and the Cauchy-Schwarz inequality yield that 

hi d v {ux-u 2 ){-Av)dS < 0{l)\\y/d^pe~£ d v {ui - u 2 )\\ L *{m + ) 

Jdn\F 



< 0(h 3 / 2 )\\e^C AmtqW ( Ul -u 2 )\\ L 2 {a) < 0(h^ 2 ) 



Here we use the boundary Carleman estimate (13. 2p and proceed similarly to 
API]) . Hence, fl4TT0|) follows. 

5. Determining the first order perturbation 
Multiplying (14 . 8 1) by h and letting h -> +0, and using (14.91) and (14. lOj) . we get 
/ (A (2) - ■ (Vp + iV^afQ^dx = 0, (5.1) 

where a^\a^ 6 C°°(f2) satisfy the transport equations (14.61) . 



Consider now (15.1 ft with <Zq = e* 1 and ag 2 " 1 = e* 2 such that 

Aw + iAii 
(Vv? + iW>) • V$ 2 + — = 0, 

a a / (5-2) 

(Vy> + iVip) • V$i + — = 0, 

and $j G C°°(Q), j = 1,2. In the coordinates (z,9), introduced in Section EJ the 
equations (15. 2p have the following form, 

2(z-z) 2(z-z) 

Hence, 

77—2 

d- z (®2 + $i) r = 0. (5.3) 

z — z 

Notice that ge® 2 with g e C^itt) such that 

{Vip + iVip) ■ Vg = in ft, (5.4) 
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also satisfies the transport equation (I4.6p . In the coordinates (z, 9), the condition 
f )5.4p reads d s g = in £7. 

Substituting = ge® 2 and = e* 1 in (15. ip . we get 

f (^(2) _ . ( Vyj + iV^ge^+^dx = 0. (5.5) 
Jn 

In the cylindrical coordinates (z,8), by (12.151) . the identity (15. 5p has the form, 

f ( A &) _ A 0)\ . ( ei + ie r )^ 2+Wl r n - 2 drdedx x = 0. 
Ja z 
Let Pg be the two-dimensional plane consisting of points (x\,r9) for 9 fixed, and 
write VLq = QnPg. We also use the complex variable z — X\ +ir, which identifies 
Pg with C. 

Taking g = g±(z) <g) g2(0) G C°°(Q), where g% is holomorphic and varying g2, we 
obtain as in [7j, for almost all 9 G § n ~ 2 , 

/ (A^ - A (1) ) • (ei + «e r )^e* 2+¥r (^ - z) n ~ 2 ^ A dz = 0, 

and hence, for all 9 G § n_2 , by continuity. Since Imz > in Q, the function g\jz 
is an arbitrary holomorphic function and therefore, we can drop the factor 1/z. 
We get 

f (A® - A (1) ) ■ ( ei + ze r )( ?1 e* 2+ * r (z - ^)"~ 2 ^ A = 0. 
By (15. 3p . we conclude that 

d- z (e^ +Wl (z - z) n - 2 ) = 0. (5.6) 

Since Imz > in ft, the holomorphic function e* 2 " 1 "* 1 ^ — z) n_2 is nowhere 
vanishing, and we can choose g\ = (e* 2+$1 (z — z) n ~ 2 )~ 1 g , where d^go = in 
fl. We get 

/ (A (2) - A {1) ) ■ { ei +ie r )g dz Adz = 0. (5.7) 

Choosing complex geometric optics solutions U2 and v as in (I4.4p and (14. 5p . where 
if) is replaced by —if), we conclude, by repeating the arguments above, that also 

/ (A {2) - A w ) ■ (ei -ie r )g dz A dz = 0, (5.8) 

where = in fl. 

Thus, since any £ G P# is a linear combination of e% and e r , choosing g = g = 1 
in (15. 7p and (I5.8p . we get 

f (A {2) - A {1) )-£dz Adz = Q for all £GP e . 
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At this point we are exactly in the same situation as the one, described in [TJ 
Section 5], see formula (5.7) there. Repeating the arguments given in that paper, 
following this formula, we obtain that d(A^ — A^) = in f2. Here and 
are viewed as 1-forms. We may notice that the arguments of Section 
5] are based upon the microlocal Helgason support theorem combined with the 
microlocal Holmgren theorem. 

Since Q is simply connected, we may write 

A (2) _ A (i) = W; ( 5 9 ) 

with * G C 5 (0). 

Our next step is to prove that $ = on O. Writing (15.71) and (15.8j) in the 
Euclidian coordinates and using (15. 9p and (I2.15p . we have 

/ (x x + ir)(Vip + «W) • (Vtygodxidr = 0, 
Jn e 

/ (xi - ir)(Vip - iVip) ■ (V*)<7 dxidr = 0. 
Thus, by (12.161) we get 

/ (d^)g dz A dz = 0, / (d z ^)g dz A dz = 0, 
J rig j \ i@ 

for all 9 G § n ~ 2 , g , g e C°°{Tr e ) such that d g g = in Q e , and <9^ = 0. By 
Sard's theorem, the boundary of Qe is smooth for almost all 9, see [3, [21] for 
more details. Thus, by Stokes' theorem, we get 



Wgodz = 0, / Vg dz = 0, 
dn g Jdn e 

for almost all 9 G S n_2 . Hence, taking go = go, we obtain that 

Vg Q dz = 0, 
and therefore, 

/ (Rety)g dz = 0, / (Im*)^^ = 0, 

Jau g 

for any g e C°°(Qe) such that <9 f go — and for almost all 9. At this point 
we are precisely in the same situation as the one, described in the beginning of 
[7J Section 6]. Repeating the arguments of that paper, we conclude that \l/ is 
constant along the connected set 8Q, and we may and shall assume that \I/ = 
along dQ. 
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Going back to (I5.ip . we get 

I (W) • (V<p + iVi))afafdx = 0, (5.10) 
Jn 

where , G C°°(Q) satisfy the transport equations (14. 6p . Integrating by 
parts in (I5.10P and using the fact that \1/ = along dQ, we obtain that 

J V ^(Atp + iAip) + (Vcp + iW) • a^a^dx = 0. 
This implies that 

J ^((Ta^a^ +a { o ) Ta ( o ) Jdx = 0, (5.11) 

where the operator T is given by (12. lip . We choose now ag G C°°(fi) being a 

solution of the equation Tag 1 '' = of the form ag = e* 1 . As for ag G C°°(fi), 
we require that Tag 2 '' = e* 2 g, where $2 G C°°(Q) is such that Te* 2 = and 
g G C°°(TL) is such that (Vy? + iVip) ■ Vg = 0. It is clear that T 2 a[, 2) = 0. The 

existence of such ag 1 ^ and a[, 2 ^ is explained in Section EJ Thus, it follows from 
@m$ that 

/ #e $2+1T adx = 0. (5.12) 
in 

In the coordinates (z,6), (15 . 12[) reads 

/ Vge** + *i(z - z) n ~ 2 dz A A dd = 0. (5.13) 
Jn 

As before, the function e* 2+<I,1 (;z — z) n ~ 2 is nowhere vanishing holomorphic in z, 
see ( I5.6p . and we shall take g = (e* 2+$1 (z — z) n ~ 2 )~ l g) g2(0), where #2 is smooth. 
Hence, (I5.13P implies that 

/ ^(x 1 ,r,6)g 2 (0)dx 1 drd9 = 0, 
Jn 

for any smooth function g2{0). We arrive exactly at the formula (6.2) of the paper 
[TJ, and arguing as in [7J, appealing, as before, to results of analytic microlocal 
analysis, we obtain that \l/ = in VL. Hence, we conclude that = in fl. 

The final step in proving Theorem 11.11 is to show that = in fi. To this 
end, substituting A^ = A^ in (14. 8 j) and letting h — > +0, we get 

/ (g (2) - q^)a^a^dx = 0- (5.14) 
Jn 
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Here we use the fact that 

_ d u (—A(ui — u-2))vdS 
an\F 



< 0(h 1/2 ) 



d v (u\ — U2)(— Av)dS 



an\F 



< 0(h 



1/21 



(5.15) 

where «2 and v, given by (I4.4p and f 14 .51) . Notice that (" 15 .151) is obtained similarly 
to (l4~9l) and (14.101) but under the assumption that = A (2) . 

We choose a^\a^ G C°°(f2) to be such that Uq = e* 1 and = e* 2 g, where 
$i, $ 2 G C°°{Q) are such that Te* 7 = 0, Te* 2 = 0, and g G C 00 ^) is such that 
(V<p + iVtp) -Vg = 0. Thus, fl5TT4"j) yields that 

( g (2) _ q W)e* 2+1 ^gdx = 0. 

Arguing in the same way as after (I5.12p . we conclude that g 1 - 1 ^ = g^ 2 - 1 . This 
completes the proof of Theorem 11.11 
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